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Abstract 

We calculate the inclusive charmless semileptonic decay width of the B meson 
using a QCD-based approach. This approach is able to account for both dy- 
namic and kinematic effects of nonperturbative QCD. The charmless semilep- 
tonic decay width is found to be enhanced by long-distance strong interactions 
by +(7± 5)% with respect to the free quark decay width. Using the resulting 
theoretical value for the charmless semileptonic decay width, an extraction of 
\V u b\ is made from the measured lifetime and inclusive charmless semileptonic 
branching ratio of b hadrons. 
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The ALEPH [|TJ and L3 || Collaborations have recently reported the first measurements 
of the inclusive charmless semileptonic branching ratio of beauty hadrons at LEP. Their mea- 
surements have been performed exploiting the different kinematic properties differentiating 
the rare b — > X u iv decay from the dominant b — > X c lv decay (the symbol I indicates either 
an electron or a muon). The inclusive charmless semileptonic branching ratio of beauty 
hadrons is measured to be 

Br(b -> X u tv) = (1.73 ± 0.55 ± 0.55) x 10~ 3 (ALEPH), 

(3.3 ± 1.0 ± 1.7) x 10 -3 (L3), (1) 

where in each collaboration the first error is statistical and the second systematic. 

The measured inclusive charmless semileptonic branching ratio combined with the mea- 
sured beauty hadron lifetime can be used to determine the Cabibbo-Kobayashi-Maskawa 
matrix element \V u b\, provided the corresponding inclusive charmless semileptonic decay 
width except for \V u b\ is theoretically predicted. This determination of \V u b\ would have a 
smaller theoretical error than a determination of it from a partial charmless semileptonic 
decay width, since, generally speaking, the latter would be more sensitive to the nonper- 
turbative QCD effects. The purpose of this work is to calculate the inclusive charmless 
semileptonic decay width of the B meson and estimate the theoretical uncertainty in this 
calculation, and then using the result to extract \V u b\. 

The inclusive charmless semileptonic decay of the B meson is induced by the underlying 
b — > u weak transition. However, experimental data are collected from decay processes at the 
hadron level, which is connected by the nonperturbative QCD interaction to the underlying 
quark-level transition. Therefore, in order to calculate the semileptonic decay width, the 
nonperturbative QCD effect has to be accounted for. A theoretical treatment of the non- 
perturbative QCD effect has been developed on the basis of the light-cone expansion 
and the heavy quark effective theory (HQET). This approach is from first principles and the 
nonperturbative QCD effect can be computed in a systematic way, leading to conceptual 
and practical improvements over the phenomenological models of inclusive decays such as 
the DIS-like parton model and the ACCMM model ||. Consequently, model independent 
predictions and a control of theoretical uncertainties become possible. 

In this approach the nonperturbative QCD effects in inclusive semileptonic B decays 
are ascribed to the distribution function /(£), characterizing the matrix element of the 
light-cone bilocal b quark operator between the B meson states. The distribution function 
/(£) can be interpreted as the probability of finding a b quark with momentum £P inside 
the B meson with momentum P. Nonperturbative QCD contributions to the inclusive 
semileptonic decay of the B meson consist of the dynamic and kinematic components. A 
crucial observation which emerges from this approach is that both dynamic and kinematic 
effects of nonperturbative QCD must be taken into account The latter results in the 
extension of phase space from the quark level to the hadron level, shown in Fig. [L] for 
the b — ► u decay, and turns out to constitute a large part of the nonperturbative QCD 
contribution to the decay width, which the heavy quark expansion approach || fails to take 
into account. 

The inclusive charmless semileptonic decay width of the B meson can be expressed as a 
convolution of the soft nonperturbative distribution function with a hard perturbative decay 
rate Yb^uiv'- 
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FIG. 1. Phase space for the b — > u inclusive semileptonic decay. The interior of the solid curve 
is the hadron level phase space (the changeableness of the mass of the final hadronic state is not 
shown explicitly). The interior of the dashed curve is the quark level phase space. 
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where M B denotes the B meson mass. We have neglected the masses of leptons and the 
u quark. The interplay between nonperturbative and perturbative QCD effects has been 
accounted for since the separation of perturbative and nonperturbative effects cannot be 
done in a clear-cut way. Thus the decay width is expressed in terms of the B meson mass 
rather than the h quark mass my,. Equation (0) includes the leading twist contribution; 
higher-twist contributions are expected to be suppressed by powers of M B (or by Sudakov- 
exponentials). 

In the /(£) — > 8(£ — m h /M B ) limit, inclusive decay rates are given by free quark decay. 
The first term in Eq. (0) corresponds to the Born-term for free quark decay. The remaining 
terms in Eq. (|2|) contain the perturbative QCD corrections, present formally as a series in 
powers of a s ji\. The perturbative QCD corrections in Eq. @ are known to first order in a s 
with the short- distance coefficient JT0[ 
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Only the parts of the higher-order perturbative corrections originated from the running of 
the strong coupling, the so-called Brodsky-Lepage-Mackenzie (BLM) corrections [|ll|], were 

The BLM corrections are presumably 



calculated to order a 2 s 



12 and to all orders O 
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dominant and may provide an excellent approximation to the full higher-order corrections. 
In our calculation, we include the second-order BLM correction with the coefficient [12|] 

cf L M = -3.22/3o , (4) 

where (3q = 11 — 2nf/3 is the first coefficient of the QCD /3-function. The second-order 
BLM correction appears to be rather large. However, this does not mean that the overall 
a 2 correction to the charmless semileptonic decay width is sizeable and a breakdown of 
perturbation theory. It was shown JI3| , P^4] | that the large value of the second-order BLM 
correction is solely due to the use of the pole quark mass in the perturbative calculation. 
The perturbation series is better behaved when the B decay rate is written by replacing the 
b quark pole mass with the T mass [|15f| . 

In order to calculate the charmless semileptonic decay width of the B meson, we need to 
know the distribution function /(£)• Due to current conservation, it is exactly normalized 
to unity with a support < £ < 1. Two additional sum rules for the distribution function 
can be obtained by using the operator product expansion and the HQET method [BJ. These 
two sum rules determine the mean value \i and the variance a 2 of the distribution function. 
The mean value is the location of the "center of mass" of the distribution function and 
the variance is a measure of the square of its width, which specify the basic form of the 
distribution function. The sum rules read H-f|: 
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A 2 = j^j- < B\hg s G aP a^h\B >, (8) 

are the HQET parameters, which parametrize the nonperturbative QCD effects on a variety 
of phenomena. The parameter A2 can be extracted from hadron spectroscopy: 

A 2 = - M 2 B ) = 0.12 GeV 2 . (9) 

The parameter Ai suffers from large uncertainty. Fortunately, from the numerical analysis 
to be discussed below we find that the result for the charmless semileptonic decay width of 
the B meson is insensitive to the variation of Ai. For our calculation we adopt the value of 



it from QCD sum rules [16 



Ai = -(0.5 ±0.2) GeV 2 . (10) 
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These properties constrain the distribution function considerably. 

Nonperturbative QCD methods such as lattice simulation could help determine further 
the form of the distribution function. The distribution function could also be extracted 
directly from experiment [[L7j]. Since these are as yet not done, we perform the calculations 
using the parametrization Q of the distribution function 

/«>= V-^+V K)g(i ~ a (n) 

where a, b, a, and (3 are four parameters, which are constrained by the sum rules (|5|) and 
(§), and N is the normalization constant. 

The distribution function /(£) seems to us not unrelated to the shape function introduced 
in |18] by resumming the heavy quark expansion. It would be interesting to discuss the 
relation of our approach to theirs. Their expressions for the decay rate are very similar, 
although somewhat different, to ours obtained in the framework of light-cone expansion. 
We can reproduce the formulas of ]Tj| by assuming £ = m^/M^. In that case, the scaling 
variable would be fixed to the definite value, instead of distributing in its entire range from 
to 1 with the probability /(£). The predictions of these theoretical approaches are subjected 
to experimental tests. 

Including both the nonperturbative and perturbative QCD contributions in the coherent 
way as described above, we are able to calculate the inclusive charmless semileptonic decay 
width of the B meson using Eq. (fj). The decay width is evaluated and the theoretical 
uncertainties are estimated as follows. 

• We study the variation of the decay width with respect to the mean value and the 
variance of the distribution function setting a = (3 = 1 in Eq. ([0]). Actually, this 
amounts to the study of the decay width as functions of m;, and Ai, since, essentially, 
the mean value of the distribution function is determined by the b quark mass and 
its variance is determined by Ai according to the sum rules in Eqs. (|5]) and (|6|). The 
results are shown in Fig. varing the b quark pole mass in the range 

m b = 4.9 ±0.15 GeV, (12) 



and Ai in the range of flTOl). The variation of m& leads to an uncertainty of 15% in 
the decay width if other parameters are kept fixed. A small uncertainty of 4% in the 
decay width results from the variation of Ai. In other words, the charmless decay 
width displays a strong dependence on the mean value of the distribution function of 
the b quark inside the B meson, but is insensitive to the variance of the distribution 
function. 

• We examine the further sensitivity of the decay width to the form of the distribution 
function when keeping the mean value and variance of it fixed, by varying the values 
of the two additional parameters a and (3 in the parametrization flTT|). We find that 
the variation of the decay width is typically at the level less than 1% if the form of the 
distribution function is changed but with the same mean value and variance, hence 
such a change has a negligible impact on the theoretical uncertainty. 
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FIG. 2. Dependence of the charmless semileptonic decay width on the theoretical input pa- 
rameters mb, Ai, and \i r . The solid (dashed) curves are for the renormalization scale \i T = rrib/2 
(fi r = 2m&). The curves with solid dots, boxes, triangles correspond to m& = 4.75,4.9,5.05 GeV, 
respectively. 



• We estimate the uncertainty due to the truncation of the perturbative series in Eq. (fj) 
by varying the renormalization scale between rrib/2 and 2r?v We also show in Fig. |2| 
the variation of the decay width due to the change in the renormalization scale. We 
observe that an uncertainty of 15% in the decay width stems from the renormalization 
scale dependence, in agreement with the estimate obtained in [13] based on a different 
method. 

This analysis implies that at present the theoretical error has two main sources: the value 
of mft (or equivalently, the mean value of the distribution function) and the renormalization 
scale dependence. The functional form of the distribution function, which is not completely 
known, is likely to cause some model dependence. However, the inclusive charmless semilep- 
tonic decay width of the B meson calculated in this approach is nearly model-independent 
since it is only sensitive to the mean value of the distribution function, which is known 
from HQET, but insensitive to the detailed form of the distribution function. It is not 
surprising that inclusive enough quantities like the total decay width are less sensitive to 
nonperturbative QCD effects. Finally, adding all the uncertainties in quadrature we find 



Y{B -> XJv) = (76 ± 16)|K 6 | 2 PS" 



(13) 



To see the impact of the nonperturbative QCD effect on the decay width, in Fig. [3] we 
compare the decay widths calculated in our approach, the free quark decay model, and the 
heavy quark expansion approach. The result in our approach shows that the nonperturbative 
QCD contributions enhance the decay width by +(7 ± 5)% with respect to the free quark 
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FIG. 3. Charmless semileptonic decay width as a function of Ai calculated in our approach 
(solid curve), the free quark decay model (dotted curve) and the heavy quark expansion approach 
(dashed curve). We take m& = 4.9 GeV and a s = 0.2. 

decay width, in contrast to the result obtained using the heavy quark expansion approach 
H (see also 0]) where a reduction of the free quark decay width by — (5 ± 2)% is found. 
The enhancement is mainly due to the extension of phase space from the quark level to the 
hadron level, indicating the importance of including the kinematic nonperturbative QCD 
effect. This result is in accord with the finding of a similar calculation of the inclusive b — ► c 
semileptonic decay width in [|J. 

The resulting theoretical value for the charmless semileptonic decay width of the B meson 
in Eq. (|i~3|) can be used to determine \V u b\. Assuming that the average charmless semileptonic 
decay width of b hadrons is the same as that of the B meson, using the average b hadron 
lifetime = 1.564 ± 0.014 ps and the inclusive charmless semileptonic branching ratio 
of b hadrons given in (|l|) measured by the ALEPH and L3 Collaborations, we extract 



\V ub \ = (3.81 ± 0.86 ± 0.40) x 10~ 3 (ALEPH), 
(5.27 ± 1.60 ±0.55) x 10" 3 (L3), 



(14) 



where in each case the first error is experimental and the second theoretical. 

In conclusion, we have extracted \V u b\ from the inclusive charmless semileptonic branch- 
ing ratio of b hadrons measured recently by the ALEPH and L3 Collaborations. The inclusive 
charmless semileptonic decay width of the B meson has been calculated in a QCD-based 
approach with theoretical uncertainties under control. This approach is able to include both 
dynamic and kinematic components of the nonperturbative QCD effects. We have shown 
that including the latter is crucial, such that the charmless decay width is enhanced, which 
leads to a decrease in the extracted value for \V u b\. It is encouraging that the nearly model- 
independent determination of \V u b\ from the inclusive charmless semileptonic branching ratio 
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of b hadrons is consistent with those determined by measuring the lepton energy spectrum 
above the endpoint of the inclusive b — > civ spectrum or by measuring the branching ratios 



for the exclusive decays B — > rriu and B — > ptv |19| . 

Currently, the error bar associated with \y u \\ determined from the inclusive charmless 
semileptonic branching ratio of b hadrons is not predominantly theoretical in origin. Both 
experimental and theoretical efforts are needed to reduce the error. We observe that the 
uncertainty from perturbative QCD is already one of the major theoretical limiting factors, 
comparable in size to another one, i.e., the uncertainty from nonperturbative QCD. To 
reduce the error due to the perturbative QCD correction, a complete calculation to order 
ot 2 s is necessary. Moreover, to reduce the theoretical uncertainty from nonperturbative QCD 
requires more accurate knowledge of the mean value of the distribution function (or the b 
quark mass). Higher- twist contributions, which are not included in our calculation, also 
deserve investigation, which may still yield numerically sizable effects of order Aqcd/^b ~ 
5%. 

In order to extract \ V u t,\, the overwhelming "background" of b — > c transitions necessitates 
the introduction of experimental cuts if statistical errors are to be minimized. The hadronic 
invariant mass spectrum appears to be a promising observable for discriminating b —>■ u signal 
from b — > c background PU| . We believe that in the future a measurement of the differential 
decay rate dT(B — > X u iv)jd^ u will be very useful for a high precision \V U \\ determination 
with a minimum overall (theoretical and experimental) error. It has been shown [17] that 
this differential decay rate is proportional to the nonperturbative distribution function. The 
measurement of the decay distribution with respect to the scaling variable £ u is of particular 



significance in a range of issues. They include [ 171] : 

• The B — > X u lv decay can be differentiated from the B — > X c iv decay even more 
efficiently than the cut on the hadronic invariant mass. 

• A model- independent determination of |V^&| can be obtained by avoiding hadronic 
complications. 

• Improved measurement of the inclusive charmless semileptonic branching ratio of b 
hadrons can be achieved. 

• The nonperturbative distribution function can be extracted directly. 
It remains to be seen whether such a measurement is feasible. 
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